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Many-body phenomena, a key interest in the investigation of bulk sohd state systems, 
are studied here in the context of the x-ray edge problem for mesoscopic systems. We 
investigate the many-body effects associated with the sudden perturbation following the x-ray 
excitation of a core electron into the conduction band. For small systems with dimensions at 
the nanoscale we find considerable deviations from the well-understood metallic case where 
Anderson orthogonality catastrophe and the Mahan-Nozieres-DeDominicis response cause 
characteristic deviations of the photoabsorption cross section from the naive expectation. 
Whereas the if-edge is typically rounded in metallic systems, we find a slightly peaked 
if-edge in generic mesoscopic systems with chaotic-coherent electron dynamics. Thus the 
behavior of the photoabsorption cross section at threshold depends on the system size and 
is different for the metallic and the mesoscopic case. 

§1. Introduction 

Many-body phenomena such as the Kondo effect or Fermi edge singularities 
(FES) have been a key interest in condensed matter physics for many years. ""^^ Moti- 
vated by the experimental progress in the field of mesoscopic physics and quantum 
chaos, especially the growing interest in many-body effects in those systems, 
we report here theoretical results on the mesoscopic x-ray edge problem. We are 
in particular interested in phenomena associated with a sudden perturbation of a 
mesoscopic system such as a quantum dot or a metallic nanoparticle. We predict 
substantial differences to the metallic case that are falsifiable in state-of-the-art ex- 
periments. 

In the x-ray edge problem,^) a sudden, localized perturbation is caused by an 
x-ray exciting a core electron into the conduction band, leaving a core hole behind. 
The response of the conduction electrons to the resulting attractive potential leads to 
Anderson orthogonality catastrophe (AOC)^^ - the overlap between the unperturbed 
and perturbed many-body wavefunctions vanishes in the thermodynamic limit. AOC 
competes with a second many-body effect known as Mahan's exciton^^ or Mahan- 
Nozieres-DeDominicis^-' (MND) response. In the metallic case, this leads to Fermi 
edge singularities, i.e., deviations from the naively expected photoabsorption cross 
section in the form of a peaked or rounded edge. More precisely, the behaviour at 
threshold is known to follow a power law, with the exponent determined by the partial 
wave phase shifts 6i at the Fermi energy in response to the sudden perturbation Vc 
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for orbital channel i^-* {cu is the x-ray energy, cuth indicates the threshold energy), 

A(a;) cx (a; - wth) ^^L'^J . (l-l) 

The two terms in the exponent have opposite signs and correspond to the MND 
response (with Iq being the optically excited channel) and to AOC, respectively. 

The many-body enhancement depends, via the dipole selection rules, on the 
symmetry relation between the core and local conduction electron wavefunction. The 
MND response will be non-vanishing only if the dipole selection rules are fulfilled. 
Assuming the local part of the conduction electron wavefunction to be of s-type, 
we distinguish between core electrons with s-symmetry (X-shell, Iq = 1) and p- 
symmetry (L2,3-shell, Iq = 0)^^ and refer to the photoabsorption threshold as K- 
or L-edge, respectively. In metals, the phase shifts are such that the if -edge is 
typically rounded whereas the L-edge is peaked. In the following, we will apply the 
usual model of a spherically symmetric potential V^^-* such that Si = for I > 0. The 
origin of the form of the FES typically observed in metals becomes then immediately 
apparent.^)' ^) 

For coherent systems with chaotic or regular dynamics we use a random matrix 
modcl*^^ or exact solutions of the Schrodinger equation, respectively, to compute the 
AOC overlap and the photoabsorption cross section. For the latter we use the Fermi 
golden rule approach introduced by Tanabe and Ohtaka.^^ Our model applies to 
nanoparticles and quantum dots with chaotic or regular (for example circular) shape, 
respectively. Comparing our results with the well-understood metallic problem, we 
find substantial changes:^)' ^'^^ (1) the finite number of particles causes AOC to be 
incomplete, (2) the sample-to-sample fluctuations of the discrete energy levels pro- 
duce a distribution of AOC overlaps, and (3) most importantly, the dipole matrix 
elements connecting the core and conduction electrons are substantially modified. 
One of our key results is that a photoabsorption cross section showing a rounded 
edge in metals will change into a slightly peaked edge on average as the size of a 
chaotic system is reduced to the mesoscopic-coherent scale. This peak is a direct 
signature of a coherent-chaotic dynamics of the conduction electrons reached in the 
mesoscopic regime: It is this property that leads to a non-vanishing dipole matrix 
element^^ and therefore to a situation that is reminescent of the L-edge behaviour. 
We will come back to this in more detail below. 

The outline of the paper is as follows. First, we consider AOC for mesoscopic 
systems subject to a sudden, rank-one perturbation. Then, we present results on the 
photoabsorption cross section in those systems. Our conclusion includes a discussion 
of possible experimental setups that allow for a verification of our prediction. 

§2. Anderson orthogonality catastrophe in mesoscopic systems 

We describe the non-interacting system by a Hamiltonian = Si o- ^*'^lo-'^«-o-' 
where ct^ creates a particle with spin a = ± in the orbital (pi{r) {i = 0, . . . , N — 1). 
The unperturbed energy levels are system specific and provide together with the 
<Pi a unique characterisation of the system. As reference point, we define the bulklike 
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system where the energy levels {ej} are spaced equidistantly and the wavefunctions 
are constant. 

We furthermore assume that the perturbing potential is a contact potential 
Vc = Vfc|rc)(rc|, with the location of the core hole and V the volume in which 
the electrons are confined. The diagonal form of the perturbed Hamiltonian is H = 
Hq + Vc = Xic] ^Ci^a, where c] ^ creates a particle in the perturbed orbital ipi{r). 
For relations between the {ej}, {Aj} and {^^i}, {i^i} we refer the reader to Refs.'^^' ^'^^ 

A remarkable property of a rank-one perturbation such as a contact potential 
is that all the quantities of interest for the x-ray edge problem can be expressed 
in terms of the {e} and {A} (or otherwise, for example the wavefunction derivative 
needed for the dipole matrix element, can be taken as independent random variables 
following a known, often a Porter-Thomas, distribution). Ignoring for now the spin 
variable, the overlap between the many body ground states with M particles of Hq 
and H, |^o) and I'P'o), can be expressed as^-' 

M-l N-1 , , , 

i=0 j=M ^ ^ J ' 

For the Fermi energy in the middle of the conduction band, the phase shift 
(5o, the perturbation strength Vc and the mean level spacing d are related through 
(5o = arctan(7ruc/c^)^^ ((^o is negative since the core potential is attractive). In our 
case, it turns out to be necessary to take boundary effects into account,^*'-* which, in 
addition to the formation of a bound state that is discussed in more detail in Ref.,^'^) 
modify the phase shifts away from the band center. This can be included simply by 
introducing a variable Vi given by 

1 1 1 iV + 0.5-i 
— = h - In . 

Vi Vc d i + 0.5 

for i € (1, A^/2), and the analogous form for i E {N/2,N). This gives rise to level- 
dependent phase shifts 6i. It is known from the metallic case that the phase shift 
at the Fermi energy determines the FES. We will now address the question to what 
extent this statement holds in the mesoscopic case. 

In Fig. [H we first discuss AOC for two mesoscopic systems of different size (with 
space for A^=100 and A^=1000 electrons in the conduction band, cf. left and right 
panels, respectively). The AOC overlap as a function of filling of the conduction 
band and the perturbation strength is shown in color scale (for reasons of symmetry, 
we also consider positive Vc in addition to the attractive Vc < describing the effect 
of the core hole). In the top row, the situation in the bulklike case is shown. Clearly, 
for otherwise equal parameters, the overlap becomes smaller for larger systems (i.e. 
closer to the thermodynamic limit). There is, however, a considerable amount of 
structure visible beyond this. The somewhat counterintuitive increase of the overlap 
with increasing filling is a property of the rank-one model that we use, in particular 
of the level-dependent phase sifts discussed above. More precisely, within this model 
the phase shift is, for Vc < 0, larger for smaller fillings, see Ref.^'^^ for details. In 
other words, the experimentally relevant phase shift at the Fermi energy depends on 
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Fig. 1. Anderson orthogonality catastrophe in bulklike systems (top row) and individual mesoscopic 
systems of diflterent sizes. The larger (smaller) system with space for up to 1000 (100) electrons 
in the conduction band is shown on the right (left). The mesoscopic system considered here is a 
circular quantum dot of radius a. The radial position rc of the perturbation is slightly different 
in the central and lower row, respectively, giving rise to different behavior of the individual 
systems. 



both the fining and the perturbation strength Vc, and one and the same phase shift 
can be realized using different sets of parameters, cf. the discussion in the context of 
Fig. [2] below. 

The central and lower row of Fig. [T]show snapshots of some corresponding meso- 
scopic cases. More specifically, a regular quantum dot of disk shape and with hard 
walls is considered. To this end, the Schrodinger equation is solved exactly for the 
energy levels and wavefunctions. The perturbation is placed at two different loca- 
tions Tc. Fluctuations characteristic for mesoscopic systems are clearly visible. In 
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Fig. 2. Distribution of the Anderson overlap \A\^ in chaotic mesoscopic systems. The same phase 
shift Sf ~ —tt/2 at the Fermi energy is realized using different parameter sets {N, M,Vc/d} 
for random matrices from (a) the COE and (b) the CUE ensemble. Once the values \A\^ are 
scaled by the system-size dependent bulklike overlap l/i^P, the distributions all coincide: The 
fluctuations of the overlap depend only on the phase shift 5f at the Fermi energy. 



particular, the overlap is not any more a monotonous function of filling: Changing 
the filling corresponds to changing the orbital at the Fermi energy. More important 
than a (model-specific) trend in the phase shift is now the distribution of energy 
levels around the Fermi energy and the amplitude of the wavefunctions at the posi- 
tion of the perturbation. This also explains the sensitivity of the Anderson overlap 
against changes of the location of the perturbing potential that are clearly visible 
when comparing the central and lower row of Fig. [TJ 

In order to gain a better understanding of the big mesoscopic fluctuations, we 
will now turn to the universal case of chaotic mesoscopic systems, regular systems 
such a circular quantum dots will be considered elsewhere. We use a random 
matrix model to effectively describe the (non-interacting) conduction electrons in 
the absence of the core hole, i.e., we assume the unperturbed energy levels {ej} 
to be the eigenvalues of a random matrix (belonging to the circular orthogonal or 
unitary ensemble, COE or CUE, respectively). The single particle wavefunctions 
will then show the characteristic Porter-Thomas probability distributions^^-* charac- 
terizing the spatial dependence of the wavefunction intensity. Ensembles of 10000 
individual chaotic systems are generated in a Metropolis algorithm; ^'^^ the joint prob- 
ability distribution for the {ej} and {Aj}, the basic ingredient for this method, was 
derived in Ref.^^^ The Anderson overlap (and later on the photoabsorption cross sec- 
tion) is computed for each realization from Eq. (|2-ip . Subsequently, average values, 
probability distributions, etc. are easily determined. 

That even in mesoscopic systems the phase shift at the Fermi energy is a physi- 
cally important quantity as known from the metallic case, becomes clear in Fig. [2l In 
Fig.[2]^a,b), probability distributions of the Anderson overlap are compared for meso- 
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scopic systems of various sizes possessing time reversal symmetry (COE case, on the 
left) or not (CUE case on the right). All parameter combinations {A'", M,Vc/d} yield 
a very similar phase shift Sp ~ —it/2 at the Fermi energy. Nonetheless, i.e., unlike 
the expectation based on the behavior of bulk systems would suggest, the probability 
distributions are rather different. But so are the reference bulklike values for the An- 
derson overlap indicated by the dashed lines (assignment is such that Af, increases in 
the order of the legend entries). Indeed, the probability distributions for the overlap 
scaled by the corresponding bulklike value convincingly coincide [lower panels (c) 
and (d)]: In the mesoscopic case, the overlap does depend on the system size N but 
the mesoscopic fluctuations are solely determined by the value of the phase shift. 

Eventually, we point out the finite probability for finding zero overlap, that 
moreover is rather different for the COE and CUE case, respectively. Referring the 
interested reader to the details given in Ref.,^'^-' the distinctive behavior traces back 
to the differences of the Porter-Thomas distribution for finding small values in the 
presence or absence of time-reversal symmetry, respectively. 

§3. Photoabsorption spectra: From rounded to peaked edge 

Next, we discuss the absorption spectra, thereby focussing especially on the K- 
edge. Our approach is based on Fermi's golden rule following the work by K. Ohtaka 
and Y. Tanabe^^ who showed that this method provides a comprehensive description 
of the x-ray edge problem. The photoabsorption cross section in the mesoscopic case 
is then obtained from (using units h = 1) 

A{u;) = 27rY,m\D\^o)\^S{Ef - E^-u) , (3-1) 
/ 

where the sum is taken over all perturbed final states (of energy Ef) connected to 
the unperturbed groundstate (Pq = n<T=± H^o^ 4 o-'^c|0) (of energy Eq) by the dipole 
operator D (cj creates the core electron in the empty band |0)). We are interested 
in processes involving the core hole; thus, the dipole operator can be written as 
D = const. ^jLo {^i'^ja^c + h.c^ . At K-edgc, the core electron wavefunction and 
the local part of the conduction electron wavefunction are both of s-symmetry, Wj 
is related to the derivative of the perturbed orbital -0^ in the direction e of the 
polarization of the x-ray through Wj = e • V'0j(rc). 

We first turn to the absorption cross section right at threshold, u) = cjth; and 
neglect the spin degree of freedom for the moment.^-*' The only possible final state 
is then = Hjlo^l^^' Without a perturbing potential, the only contribution is 
the direct process wmc^m^c- In the presence of a perturbation, however, the new and 
old orbitals are not identical. This implies that the so-called replacement processes, 
terms with j < M, also contribute coherently, giving^) 



\{^fo\D\^^)\^ oc \wmA\ 



^ wmA 

1=0 



(3-2) 
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where Aj is defined by generalizing Eq. (12-ip with level i (< M) replaced by M. 
Since for chaotic systems the derivative of the wavefunction, Vk~'^ x \ Vgipj\'^, is known 
to have Porter-Thomas fluctuations uncorrelated with the wavefunction itself/'^-* we 
can proceed as for the overlap to construct the distribution of \{^fo\D\(!>Q)\'^ . 

Away from threshold, part of the x-ray energy can excite additional electrons 
above the Fermi energy in so-called shake-up processes. Their contribution is a 
straightforward generalization of Eq. (j3-2p . Although the number of these processes 
grows in principle exponentially with the energy of the x-ray, only few shake-up pro- 
cesses contribute significantly to the photoabsorption. Shake-up processes involving 
more than three shake-up pairs can safely be neglected^^^ as was also found pre- 
viously.^^^ The spin of the electrons is taken into account by including the AOC 
contribution due to the additional electronic channel. 
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Fig. 3. Photoabsorption spectra for chaotic-coherent mesoscopic systems and bulkhke metals in 
comparison. Whereas the L-edge is peaked in both cases, the typicaUy rounded metaUic K-edge 
evolves into a peaked signature when the system size is reduced to the chaotic-coherent scale. 



The result for the photoabsorption cross section at the K-edge of a mesoscopic 
system is shown in Fig. [3] (open circles). For comparison, the bulklike curves at 
K- and L-edge (filled triangles and squares, respectively) are also provided. They 
are obtained assuming equidistant energy levels and constant dipole matrix elements 
that explicitly depend on the symmetry relation between the optically active channel 
and the core electron.^) As discussed above, this typically leads to a rounded K-edge 
(vanishing dipole matrix elements, only AOC contributes) and a peaked L-edge (the 
MND response, being linear in the phase shift, overcompensates AOC). Remarkably, 
we find a behavior reminiscent of such a peaked edge at the mesoscopic K-edge that 
is, on average, slightly peaked. This striking difference has its origin in the chaotic- 
coherent dynamics of the electrons in generic (ballistic) mesoscopic systems such as 
quantum dots or metallic nanoparticles. 



8 



Martina Hentschel et al. 



§4. Conclusions 

The central result of this work is that changes in the dynamics of electrons in 
a Fermi sea may imply characteristic changes in the many-body response of the 
system. They occur, e.g., as a result of a systematic reduction of the system size 
from the bulklike-metallic to the mesoscopic-coherent scale. One possibility to make 
these changes visible is through the photoabsorption cross section in response to 
the sudden creation of a localized perturbation following the excitation of a core 
electron. In particular, a typically rounded K-edge should develop into a slightly 
peaked edge when the system size is sufficiently reduced to induce chaotic-coherent 
dynamics of the electrons. This signature, marked by the arrow in Fig. [3l is an 
effect of the coherent confinement in the chaotic system where the dipole matrix 
element at K-edge is determined by the derivative of the wavefunction that, unlike 
the bulklike case, is independent from the wavefunction itself. Most importantly, it 
will take non-zero values on average, and consequently lead to a signature that is 
qualitatively comparable to the metallic L-edge behavior. 

Although the effect of the transition to a rounded edge seems to be rather small 
and requires resolution of the x-ray energy on the order of the mean level spacing 
(cf. Fig. [3|), such an x-ray absorption experiment using metallic nanoparticles should 
become possible in the near future. Using nowadays technology, we suggest exper- 
iments using arrays of quantum dots. The excitation would then not occur by an 
x-ray and from the core level, but rather by radiation from a microwave laser and 
from an impurity level specifically introduced by doping in between the valence and 
conduction band: The physics that we describe here, namely the sudden perturbation 
of a Fermi sea of electrons by a localized potential, is the very same. The available 
energy resolution and manageability allow one, in principle, to determine the average 
values of the photoabsorption cross section and the signature of a coherence peak at 
the -ftT-edge threshold. 
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